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I We present a theoretical and analytical curve with reproduce essential features 

of the frequency distributions vs. diameter, of the 42,000 crater contained in the 
Barlow Mars Catalog. The model is derived using reasonable simple assumptions 
i/^ [ that allow us to relate the present craters population with the craters population 

^ I at each particular epoch. The model takes into consideration the reduction of the 

O ! number of craters as a function of time caused by their erosion and obliteration, 

and this provides a simple and natural explanation for the presence of different 
slopes in the empirical log- log plot of number of craters (N) vs. diameter (D). 
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c/3 . 1. Introduction 

The present impact crater size frequency distribution, N is the result, on one hand, of 
^ I a rate of crater formation, 0, and, on the other hand, the elimination of craters, as time 
I goes by, due to effects like erosion and obliteration. Therefore if we want to understand 
the crater formation history we will need to know how these forming and erasing factors 
combine to create A^. Thus, in this work the above problem is analyzed, and in section 2 we 
find that N can be expressed in terms of (p and the fractional reduction of craters per unit 
of time, C. Then, a simple model is discussed that describe the crater size distribution in 
Mars data, collected by Barlow (Barlow 1988), where it is assumed that is independent of 
time. The above model is realistic, since according to several investigations has remained 
nearly constant for the last 3 to 3.5 billion years (Hartmann 1966b; Neukum 2001, 1983; 
Ryder 1990). The simplest interpretation of this model implies that and C are given as 
the following inverse power of the diameter, D, of the crater: oc C oc In section 
3 the model is applied to craters data on Earths, and it is concluded that also in our planet 
C oc This result is interpreted to mean that on Mars and Earth we have C ~ volume' 
or equivalently the crater mean life = ^ oc Volume oc h, with h oc Investigations 



- 2 - 



of geometric properties of Martian impact craters reflect values of the average height h{D) 
consistent with the above conclusion. 



In what follows we will present theoretical and analytical curves which will reproduce the 
essential features of the martian crater-size frequency distribution empirical curves (Figure 
1), based on Barlow's (1987) database of about 42,000 impact craters. The models will be 
derived using reasonable simple assumptions, that will allow us to relate the present crater 
population with the crater population at each particular epoch. 

To this end, let AN(D, r) represents the number of craters of diameter D ± ^ formed 
during the epoch r ± where we are assuming that AD and Ar are sufficiently large that 
is justified treating AA^ as a statistical continuous function, but, on the other hand, they 
should be sufficiently small (^ <^ ^ 1) to be able to treat them as differentials in the 

following discussion. This initial population will change as time goes on due to climatic and 
geological erosion, and the obliteration of old craters by the formation of new ones. Then, 
we expect that the change in AN during a time interval dt will be proportional to itself and 
dt: 



where C is the factor that takes into account the depletion of the craters, and should be 
a function of the diameter, since the smaller a crater is the most likely it will disappear. 
Furthermore, C could also depend on time however, we will ignore such changes here, which 
we believe is a good starting approximation to the general problem. It is easy to integrate 
Equation (1) in time to obtain: 



2. 



Theoretical Models for the Observed Data 



d{AN) = -CANdt, 



(1) 



AN{D,t) = AiV(D,y Exp [-Ct„], 



(2) 




(3) 



Equation (2) gives the number of craters, as a function of D, observed at time t, that were 
produced at the time interval t„ ± Therefore the total contribution to the present (t=0) 
population due to all the epochs t„ is: 
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N{D) = AiV(L',t„) Exp [-Cr„], 



(4) 



n 



or in the continuous limit Ar^O, 



N{D) = / (j){D,T)Exp[-CT] dr, 



(5) 







where 



4>{D,t) = hniAr^o 



AN{D,Tn) 



(6) 



At 



4>{D, t) is the rate of crater formation of diameter D at the epoch r, and r/ is the total time 
of crater formation. 

In the next section we will determine the function C{D) and for a model where we 
assumed that the rate of crater formation, 0, is independent of r. 



Investigations of the time dependance of cratering rate of meteorites have concluded 
(Hartmann 1966b; Neukum 2001, 1983; Ryder 1990) that the impact rate went through a 
heavy bombardment era that decayed exponentially until about 3 to 3.5 Gy, and since then 
has remained nearly constant until the present. Therefore, for surfaces that are younger than 
3 to 3.5 Gy we can reasonably assume that (j) is independent of r, and hence from Equation 
(5) immediately obtain 



We then find that the simplest model that essentially reproduces the data in Figure 1, for 
D > 6 km, is given by Equations (8)and (9): 



3. (f) Independent of r 



N{D) 



[l-Exp{-CTf}]. 



(7) 



C{D) 



<f>{D) 



3.55 X 10^ 




-4- 



O 
O 

o 




1 2 



20 50 200 



D (Km) 



Fig. 1. — Log-log plot of number of craters (N) vs. diameter (D) in Mars. Following (Ncukum 
2001) the number of craters per kilometer squared were calculated for craters in the diameter 
Dl < D < Dr, where and Dr are the left and right bin boundary and the standard bin 
width is Dr/Dl = 2^/2. 
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(9) 



We see that the theoretical curve (7), shown in Figure (2), differs significantly from the 

observed curves for D less than about 6 km. However, according to Barlow (Barlow 1988) 
the empirical data is undercounting the actual crater population for D less than 8 km, and 
therefore no meaningful comparison is then possible between models and data for this region 
of small craters. 

Equation (2) implies that the fraction of craters of diameter D formed at each epoch r 
that still survive at the present time r = is given by: 



AAr(L>,0) 



— Exp [— Cr] fti Exp 



57 

'd 



2.5 



(10) 



and thus we have that the mean life for craters of diameter D, Tmean, is 



-mean ^ ' \ 57 



D 



2.5 



Hence, craters with D ^ 57 km have r„ 



Tf, while 



Tmean » Tf , D » 57 km, 
Tmean « T/ , D « 57 km. 

The region D » 57 km is approximately described by the limit of Equation (7) when 
L> ^ oo: 

3.55x10^ 

j.m^iV = 0r, = ^^, (13) 

which corresponds to a straight line of slope -4.3 in a Log N vs Log D plot, and that would 
be the form of Equation (7) in the absence of erosion and obliterations (C~0). Hence,wc 
have that the bending of the empirical curve (Figure 1) for D < 57 km is explained in this 
model as the result of the ehmination of smaller craters as they get older. We also see from 
Equations(13) that when the effect of C can be ignored we have N — (j) (D) Tf, and therefore 
the actual crater density N is proportional to the age of the underlying surface Tf. On the 
other hand, when for smaller craters Exp [— C Tf] « 1 we will have from (7) that 
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Fig. 2. — Comparison of theoretical model with the empirical log-log plot of number of 
craters (N) vs. diameter (D) in Mars. 
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N{D) ^ = 0(L')Wn, (14) 

and in this limit the crater density N{D) is proportional to the survival mean life, Tmeaiu 
of the craters of size D. Thus, when saturation occurs and hence N is independent of 
Tf, we have, instead, that N is proportional to Tmcan- This feature is called by Hartmann 
(Hartmann 2002) "Crater retention age" , and in Mars this effect shows, according to this 
model, in craters smaller than about 57 km. 



4. Application to Earth 

The model given by Equations (7), (8), and (9) assumed a simple polynomial form for 
and C, however, alternative models can be also considered. For instance, by assuming that 



1.43 X 10^ 

N = 0T/ = 



^1.8 



1 - Exp 



-2.48 X 10^ 



(15) 



we will reproduce the Mars crater data, exactly as in model given by Equations (7), (8), (9) 
but now with C — , and the change in slope in Figure (1) around D ^ 57 km will now 
be interpreted as intrinsic behavior of (p{D) rather than due to the erosion and obliteration 
of smaller craters. How can we then discriminate between these two alternative views?. We 
sec that in the model given by Equation (7) the fraction of craters of a given diameter, D, 
produced at a time r, decreases with time according to Equation (10) as 

while in the model of Equation (15) this fraction is independent of time. Therefore we can 
put to test the validity of Equation (16) by studying crater size frequency distributions as a 
function of time. This is possible to do in our planet, and in this section we will investigate 
the consistency of the hypothesis (16) with the Earth craters data. 

Thus consider the average diameter of craters observed today that were formed during 
a given time r ± ^, which is given, according to Equation (16), by 

D = ^ - 17) 
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Assuming that C and behave in the form 

A 



, A — const, m — const, (18) 



C = , B — const, I — const, (19) 
we can rewrite Equation (17) in the form (Appendix) 

D = BT aTT, (20) 

where 



r 



' 1)1—2 



a = ' ) (21) 

r(^) 



and 



oo 

r(n) = / U^'-^e-^dU (22) 
Jo 

is the Gamma function. Equation (20) can be rewritten as 

Log!) = y Logr + LogBT a, (23) 

t 

which represents a hnear relation between Log 5 and Logr with slope y. In Figure (3) 
we plot LogD vs Logr from data of crater size vs r on Earth, and the straight line best 
fitting gives I — 2.5, which is the value determined for model (7) for Mars. This result is 
interpreted as follows. If we assume that, as expected, Tmean is a function of the volume of 
the crater, V, that decreases with decreasing V, then it is reasonable to expand it in terms 
of powers of V, and thus we will have 

r = i = aiF + 02^' + a^V^ + . (24) 



Furthermore, for sufficiently small volumes we would have, as a good approximation to C, 
that 
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Fig. 3. — Average diameter (D) vs. average age (f) for terrestrial craters. The bin size 
increases as 2?. The slope of the straight line best fit (0.40) correspond to I — 2.5. 
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l-^aiV ^ aiD'^h, (25) 



where we are writing 



V ^ D'^h, (26) 

with h as the average height of the crater of size D. The comparison of Equation (25) with 
Equations (19), with / = 2.5, imply that 



h - ConstD^ (27) 

which is a prediction that can be investigated, and we have found that indeed Equation (27) 
is consistent with results from studies of impact crater geometric properties on the surface 
of Mars, by J.B. Garrin (Garrin 2002). 

Therefore it appears that the age distribution of craters on Earth favor the simple model 
considered for Mars, where there is an erosion and obliteration factor C with the approximate 
form 



^ Const , , 

C (28) 

It is also suggested here that the above behavior for C follows from a relation of the form 



_ Const _ Const, 

V ~ D^h' ^ ^ 



with 



h (X D2 (30) 

Further investigations and observations of the crater data on the terrestrial planets, the 
moon and the asteroids are necessary for additional tests of the validity of the model (7) and 
its interpretation. 



5. Appendix 



Lets define 



or, equivalently 



Then we have 



u 



and therefore Equation (17) becomes 



D = {Bt) 



1 ru'^-'e-^dU 



j^U-^-^e-udU 
where a is the ratio of the following gamma functions: 
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